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The aim of the study was to implement a mathematical model to simulate two-phase anaerobic digestion 
(TPAD) process which consisted of an anaerobic continuous stirred tank reactor (CSTR) and an upflow 
anaerobic sludge blanket (UASB) reactor in series treating traditional Chinese medicine (TCM) wastewa¬ 
ter. A model was built on the basis of Anaerobic Digestion Model No. 1 (ADM1) while considering com¬ 
plete mixing model for the CSTR, and axial direction discrete model and mixed series connection model 
for the UASB. The mathematical model was implemented with the simulation software package MAT- 
LABTM/Simulinks. System performance, in terms of COD removal, volatile fatty acids (VFA) accumulation 
and pH fluctuation, was simulated and compared with the measured values. The simulation results indi¬ 
cated that the model built was able to well predict the COD removal rate (-4.8-5.0%) and pH variation 
(-2.9-1.4%) of the UASB reactor, while failed to simulate the CSTR performance. Comparing to the mea¬ 
sured results, the simulated acetic acid concentration of the CSTR effluent was underpredicted with a 
deviation ratios of 13.8-23.2%, resulting in an underprediction of total VFA and COD concentrations 
despite good estimation of propionic acid, butyric acid and valeric acid. It is presumed that ethanol pres¬ 
ent in the raw wastewater was converted into acetic acid during the acidification process, which was not 
considered by the model. Additionally, due to the underprediction of acetic acid the pH of CSTR effluent 
was overestimated. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

Traditional Chinese medicine (TCM), a splendid culture of Chi¬ 
nese nation for thousands of years, occupies the very important 
status in the medicine domain due to its unique effectiveness. In 
recent years, TCM industry rapidly develops, which results in a cor¬ 
respondingly increase in TCM wastewater production. It is esti¬ 
mated that approximately half of the TCM wastewater was 
discarded directly without specific treatment (Shi et al., 2005). 
The composition of TCM wastewater was determined by raw mate¬ 
rials and production process, mainly coming from the washing 
water of raw material, the residual liquid of original medicine 
and floor washing water (Jiang et al., 2006; Zhu, 2007; Zhao 
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et al., 2007). TCM wastewater contains a lot of various natural or¬ 
ganic pollutants, such as carbohydrates (mainly polysaccharides), 
alkaloids, phenols,alcohols, amino acid, lignin, proteins, pigments 
and their hydrolysates, in the form of water-soluble, colloidal, sus¬ 
pended particles (Li et al., 2006; Zhu, 2007). Therefore, TCM waste- 
water is a kind of high concentration organic wastewater with 
complicated composition, which will cause the serious pollution 
to the environment if emitted directly. 

Due to the high COD concentration in pharmaceutical wastewa¬ 
ters, attempts have already been made to work with anaerobic pro¬ 
cesses (Yang et al., 2003; Goblos et al., 2008). It has been noted that 
the acid inhibition usually occurs in the one-phase anaerobic di¬ 
gester because of the differences in the rates of acidogenesis and 
methanogenesis, which depresses the methane yield and the sta¬ 
bility, resulting thereby in low treatment efficiency. The different 
growth rates and pH optima for acidogenic and methanogeic 
organisms, and thus different requirements regarding reactors, 
have led to the development of two-phase anaerobic digestion 
(TPAD) processes (Liu, 1998). Currently, research on TPAD 
processes was concentrated on the evaluation of operating condi¬ 
tions, biological characteristics, and process performance, but not 
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the model development. Model study might be favorable to com¬ 
prehensively evaluate the impacts of all process variables on the 
performance of TPAD processes. 

Many models had been developed for anaerobic digestion pro¬ 
cesses, but with increasing complexity of the advanced digestion 
technologies, more complex models that can represent the impacts 
of changing environments on chemical and microbial species were 
required (Parker, 2005). Recently, there had been a move by the 
International Water Associations (IWA) Task Group for mathemat¬ 
ical modelling of anaerobic digestion processes to develop a com¬ 
mon model (Batstone et al., 2002a). ADM1 model has been used 
as a common basis for further development or as a common 
platform for dynamic simulation of different anaerobic processes. 
Blumensaat and Keller (2005) used ADM1 model to implement 
the anaerobic two-stage digestion of sewage sludge, a thermophilic 
pre-treatment followed by a mesophilic main treatment stage. The 
accuracy of the optimized parameter sets has been assessed 
against experimental data from pilot-scale experiments. Under 
these conditions, the model predicted reasonably well the dynamic 
behaviour of a two-stage digestion process in pilot-scale. Consider¬ 
ing microbial storage, Toshio et al. (2007) demonstrated that the 
modified model adequately described anaerobic sequencing batch 
reactor ASBR dynamics in the degradation of trehalose during a 24- 
h cycle. There were also many implementations of this powerful 
tool that have been tested and proved their success in simulating 
the anaerobic digestion of several organic wastes (Fuentes and 
Nicolas, 2008; Fezzani and Ridha, 2008). 

However, application of ADM1 model to the two-phase anaero¬ 
bic process, especially for the treatment of drugs manufacturing 
with two-phase anaerobic process wastewater remained limited. 
The objective of this study was to make ADM1 model modified 
based on the process configurations and then to simulate experi¬ 
mental results of two-phase anaerobic process treating TCM 
wastewater with the modified ADM1 model. 

2. Methods 

2 A. Pilot-scale setup 

A schematic diagram of pilot-scale anaerobic process is shown 
in Fig. 1. The system consisted of two anaerobic reactors, that is, 


a continuous stirred tank reactor (CSTR) and an upflow anaerobic 
sludge blanket anaerobic reactor (UASB), working as acidogenic 
conditions and mehtanogenic conditions, respectively. The CSTR 
had a working volume of approximately 12 m 3 with an internal 
diameter of 2.2 m and a height of 3.2 m. A blender was equipped 
for culture mixing. The UASB was built cuboidally with a square 
section of 9 m 2 and a height of 6.5 m, providing a working volume 
of approximately 55 m 3 . The top half of the reactor was filled with 
wave-shaped epoxy glass cells (cell thickness, 1 mm; crest height 
of wave, 45 mm; crest distance of wave 100 mm; slant distance 
of wave, 130 mm; and specific surface area, 360 m 2 m -3 ). A 
three-phase separator was installed upside the reactor to prevent 
biomass washout. The HRTs of CSTR and UASB were maintained 
consistently at 12 h and 55 h, respectively. The reactor tempera¬ 
tures were maintained through auto-controlled heat exchangers 
to be 35 °C. 

2.2. Wastewater 

The wastewater, obtained from a local TCM company (Harbin, 
China) was used as the feedstock of the CSTR-UASB system. The 
wastewater was generated from the processes of washing, distilla¬ 
tion, separation, concentration, purification washing, refining, tor- 
refaction, and moulding. Its biochemical oxygen demand (BOD 5 ) 
was estimated to be around 4200 mg/L and COD around 17000 
mg/L, indicating a low biodegradability (BOD 5 /COD « 0.25). Vola¬ 
tile fatty acids (VFA), including acetic, propionic and butyric acids, 
and ethanol accounted for about 3.65% of total wastewater COD, 
with the averages as follows: acetic acid (241 mg/L), propionic acid 
(63 mg/L), butyric acid (68 mg/L) and ethanol (358 mg/L). 

2.3. Analytical methods 

The TCM wastewater and effluent of CSTR and UASB were 
collected daily for the analysis of COD. Measurement of COD was 
performed directly on the sampled slurry according to Standard 
Methods (APHA, 1995). VFA and ethanol were determined by gas 
chromatography (GC) as described in the previous work (Ren 
et al., 2006). The GC (GC-122, Shanghai Analytical Apparatus 
Corporation, China) was equipped with an injector (220 °C), a 
hydrogen flame ionization detector (220 °C) and a 2 m 


Elevated water tank 


Methane 



Fig. 1. Scheme of the TPAD system. 
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(length) x 5 mm (ID) stainless steel column packed with supporter 
of GDX-103 (60-80 meshes). The operation of the stainless steel 
column was amenable to a temperature programming process 
within 100-200 °C (N 2 carrier at a flow rate of 50 mL/min, and 
combustion gas H 2 at 50 mL/min, combustion-supporting gas 0 2 
at 500 mL/min). The retention times were of 0.474 min for ethanol, 
0.690 min for acetic acid, 1.412 min for propionic acid, 3.020 min 
for butyric acid and 6.727 min valeric acid. Before measurement, 
the sample was filtrated through a 0.45 pm filter and acidified with 
6 M HC1. Two pL acidified sample was injected manually for anal¬ 
ysis. ORP and pH were monitored daily by ORP determinator (ORP- 
412, Cany Precision Instruments Co., Ltd, China). 

3. Mathematical model 

3.1. Substrate degradation and main hypotheses of modeling 


cause the microorganism content is neglected in influent and 
effluent flows, the calculation of the state variables is shown in 
Eqs. (1) and (2). 


Qin^ijn Qout^i + ^ ^ Pj Vij 

j= 1-19 

(i) 

~% X i + ^2 Pj V ij 

(2) 


Where 19 pjZ4/ is the sum of the biochemical rate coefficients 
( Vij ) multiplies kinetic rate equations (p-) for process; V is the liquid 
reactor volume; q in is flow rate of influent; q out is flow rate of efflu¬ 
ent; q w is flow rate of sludge discharge; Sij n is the input concentra¬ 
tion of the soluble components; 5, is the concentration of soluble 
components and Xi is the concentration of the particulate in the 
reactor. The components are shown in Table 1. 


It is recognized that a number of the conversion processes that 
are active in anaerobic digestion can be inhibited by the accumula¬ 
tion of intermediate products such as molecular hydrogen, ammo¬ 
nia or by extremes of pH. In this study, it is assumed that all 
microbially mediated substrate conversion processes are only 
inhibited by extremes of pH. Fluid-gas exchange process is ne¬ 
glected. The microorganism content is neglected in influent and 
effluent. These simplifying assumptions are introduced to make 
the model workable, although they do not completely reflect 
reality. 

3.2. Main mathematical model equations 

In ADM1 model, cellular kinetics is described by three expres¬ 
sions: uptake, growth and decay. This model includes three overall 
biochemical (cellular) steps, this is, acidogenesis (fermentation), 
acetogenesis (anaerobic oxidation of organic acids) and methano- 
genesis. In a two phase anaerobic process, acidogenic-phase diges¬ 
tion process includes acidogenesis and acetogenesis steps, while 
methanogenic-phase process involves methanogenesis. 

Matrix form is adopted in ADM1 model to describe the bioki¬ 
netics model (Batstone et al., 2002a). An important advantage of 
matrix form is able to easily understand dynamic process of the 
various components used in the model. The matrix includes bio¬ 
chemical rate coefficients () and kinetic rate equations (p-). Be¬ 


Table 1 

Model component. 


Component 

Name 

Value 

Unit 

S su 

Monosaccharides 

Simulation 

kgCOD/m 3 

S a a 

Amino acid 

Simulation 

kgCOD/m 3 

Sfa 

Long chain fatty acid 

Simulation 

kgCOD/m 3 

Sva 

Valerate acid 

Simulation 

kgCOD/m 3 

Sbu 

Butyric acid 

Simulation 

kgCOD/m 3 

Spro 

Propionic acid 

Simulation 

kgCOD/m 3 

Sac 

Acetic acid 

Simulation 

kgCOD/m 3 

Sh2 

Dissolved H 2 

Simulation 

kgCOD/m 3 

Sch4 

Dissolved CH 4 

Simulation 

kgCOD/m 3 

Si 

Soluble inerts 

Simulation 

kgCOD/m 3 

Xc 

Composites 

Simulation 

kgCOD/m 3 

Xch 

Carbohydrates 

Simulation 

kgCOD/m 3 

X pr 

Proteins 

Simulation 

kgCOD/m 3 

x H 

Lipids 

Simulation 

kgCOD/m 3 

Xsu 

Sugar degraders 

Simulation 

kgCOD/m 3 

Xaa 

Amino acid degraders 

Simulation 

kgCOD/m 3 

Xfa 

LCFA degraders 

Simulation 

kgCOD/m 3 

X c4 

Valerate and butyrate degraders 

Simulation 

kgCOD/m 3 

Xpro 

Propionate degraders 

Simulation 

kgCOD/m 3 

Xac 

Acetate degraders 

Simulation 

kgCOD/m 3 

x h2 

Hydrogen degraders 

Simulation 

kgCOD/m 3 

Xi 

Particulate inerts 

Simulation 

kgCOD/m 3 


h=h=h = exp 


/ pH-pH ut \ 2 
V pH UL -pHj 


( 3 ) 


To reduce the complexity of the model, only pH inhibition mecha¬ 
nisms are considered, as shown in Eq. (3). 


32A. The balance equations of CSTR 

In the present study, the CSTR is considered to be completely 
mixed through continuous stirring, and the frequency and inten¬ 
sity of mixing are sufficient to prevent the formation of large sub¬ 
strate gradients in the reactor. Based on the rate equations matrix 
of ADM1, the following balance equations of CSTR can be derived, 
as shown in Eqs. (4-20) (Appendix A). 


3.2.2. The balance equations of UASB 

According to the distribution of sludge concentration along the 
reactor height, UASB reactor was divided into three internal sec¬ 
tions, i.e., bed section, blanket section and settler section. In order 
to simple the flow model, the axial direction discrete model and 
mixed series connection model were introduced here. Discrete 
model, one kind of revised ideal model, mainly describes non-ideal 
reactors, particularly applicable to the non- back-mixing systems. 
This model has assumptions as follows: (1) each section vertical 
to flow direction has homogeneous radial concentration; (2) axial 
direction discrete coefficient D L is constant, regardless of operating 
time and reactor position; (3) substrate concentration is continu¬ 
ous functions of flow position. Two equations [Eqs. (21) and (22)] 
representing inflow and outflow states can de derived. 


inflow: 


qS + D 


d 

L n 



( 21 ) 


outflow : 



D dS1 
° L dl 


( 22 ) 


Where A is the section area of reactor, l is the reactor height and V is 
the reactor volume.As high order differential items are neglected, a 
balance equation [Eq. (23)] can be derived. 


dS 

dt 




(23) 


Boundary conditions for Eq. (23) are described as follows: 

1 = 0, q S in =qS +in -D L (jfj (24) 

'-L- {tD-0 < 25 > 
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According to the mixed series connection model, the actual reactor 
is divided into several equal volume ideal-mixing sectors. The mass 
balance of those sectors is expressed by Eqs. (26) and (27). 


VdSi_ 
n dt 


<?%,+ Pj v u 

j= 1 — 19 


( 26 ) 


E Pj v ‘J ( 27 ) 

1 j=l—19 

By combining the rate equations matrix of ADM1 and the flow 
model, the balance equations of UASB are derived, as shown in 
Eqs. (28)-(38) (Appendix B). Boundary conditions for the above 
material balance equations are given in following equations [Eqs. 
(24) and (25)]: 


3.3. Model parameters estimation 


In ADM 1, there are three main parameters: stoichiometric coef¬ 
ficients, equilibrium coefficients, and kinetic parameters. Many of 
the stoichiometric coefficients, equilibrium coefficients as well as 
kinetic parameters are assumed to be fixed due to their low vari¬ 
ability and sensitivity and taken from the scientific and technical 
work reported by the IWA task group and Chynoweth et al. (Chy- 
noweth et al., 1998a; Batstone et al., 2002b; Masse et al., 1996, 
1997; Masse and Droste, 1997, 2000; Stumm and Morgan, 1996). 
The reactor parameters were taken from the experimental set¬ 
ups for both datasets and applied to the model. The operating tem¬ 
peratures were 35 °C for CSTR and UASB. The HRTs of CSTR and 
UASB were maintained consistently at 12 and 55 h, respectively. 
All simulations were run to steady-state using a constant input 
including a constant flow rate that was given through the initial 
conditions of the dynamic experiments, to avoid numerical prob- 


Table 2 

Model parameters with low sensitivity and low variability from the literature. 


Parameter 

Name 

Value 

Unit 

A. Hydrolysis parameter 

fchyd-cH Hydrocarbon hydrolysis coefficient 

8 

d- 1 

^hyd-PR 

Proteins hydrolysis coefficient 

8 

d- 1 

^hyd-Li 

Lipid hydrolysis coefficient 

8 

d- 1 

B. Kinetic parameters 

k m -s U Monod largest absorption rate of sugar 

30 

COD/(CODd) 

Ks-su 

Half saturated coefficient of sugar 

0.5 

kgCOD/m 3 

Y su 

Yield of sugar 

0.1 

COD/COD 

/<m-aa 

Monod largest absorption rate of amino 

50 

COD/(CODd) 

K s-aa 

acid 

Half saturated coefficient of amino acid 

0.3 

kgCOD/m 3 

Yaa 

Yield of amino acid 

0.08 

COD/COD 

fcm-fe 

Monod largest absorption rate of long 

6 

COD/(CODd) 

K S -fa 

chain fatty acid 

Half saturated coefficient of long chain 

0.4 

kgCOD/m 3 

Yb 

fatty acid 

Yield of long chain fatty acid 

0.06 

COD/COD 

k m -c+ 

Monod largest absorption rate of butyric 

20 

COD/(CODd) 

«s-c t 

acid 

Half saturated coefficient of butyric acid 

0.3 

kgCOD/m 3 

y c +4 

Yield of butyric acid 

0.06 

COD/COD 

Ypro 

Yield of propionic acid 

0.04 

COD/COD 

Vac 

Yield of acetic acid 

0.05 

COD/COD 

k m -H 2 

Monod largest absorption rate of H 2 

35 

COD/(CODd) 

Yh 2 

Yield of H 2 

0.06 

COD/COD 

C. Stoichiometric coefficients 

/ bu su Butyric acid generated by sugar 

0.59 


/ac, su 

Acetic acid generated by sugar 

0.45 

- 

/pro, su 

Propionic acid generated by sugar 

0.20 

- 

/ac, aa 

Acetic acid generated by amino acid 

0.6 

- 

/pro,aa 

Propionic acid generated by amino acid 

0.02 

- 

/bu, aa 

Butyric acid generated by amino acid 

0.2 

- 

/v a, aa 

Valeric acid generated by amino acid 

0.3 

- 


lems that could occur when starting a simulation with dynamic in¬ 
put (Pauss et al., 1990). The prior simulations were run over an 
integration period of 120 days to find the effective steady state. 
Carbon contents were recalculated by implementing the balance 
terms in the rate equation matrix. The mass balance has been 
checked for COD and carbon contents using the Excel spreadsheet 
developed by Batstone et al. (2002a). These specific biochemical 
process rates and coefficients are tabulated in Table 2. 

Parameters with high sensitivity and high variability are esti¬ 
mated using experimental data of Masse et al., 1996,1997; Masse 
and Droste, 2000; Siegrist et al., 2002; and the selected optimiza¬ 
tion methods (Masse et al., 1996, 1997; Masse and Droste, 1997, 
2000; Siegrist et al., 2002). The nonlinear constrained optimization 
method was implemented using MATLAB 6.5 optimization toolbox. 
The two-parameter optimization around optimum using the secant 
method is implemented using Aquasim 2.1 d version. Estimation 
procedures were applied for the following parameters: the dissoci¬ 
ation coefficient, k dis ; Monod largest absorption rate of acetic acid 
and propionic acid, k miac , k m> p ro ; Half saturated coefficient of acetic 
acid, propionic acid and hydrogen, K sac , I< Sy pro , I< Sy H2 . Estimation re¬ 
sults of the six parameters were shown in 4.1 sections. 

3.4. Ordinary differential equation solutions 

Anaerobic digestion processes was formed the work of which 
led to the Anaerobic Digestion Model No. 1 (Batstone et al., 
2002b). ADM1 is a highly complex model, characterized by 19 bio¬ 
chemical conversion processes and 24 dynamic state variables. For 
this paper calculations were executed with the software MATLAB 
version 6.5 (release 13) and Simulink 5. 

MATLAB/Simulink system provides ordinary dynamic equation 
solver systems (Copp et al., 2002; Rosen and Jeppson, 2002. Solver 
ode23, a one-step solver, is an implementation of an explicit Run- 
ge-Kutta (2,3) pair of Bogacki and Shampine. It could be more effi¬ 
cient than ode45 in cases of crude tolerances and moderate 
stiffness. Solver ode45, a one-step solver, is based on an explicit 
Runge-Kutta (4,5) formula, the Dormand-Prince pair. It could be 
considered as the best function to apply as a “first try” for most 
problems. Solver odel5s, a multi-step solver, is a variable order 
solver based on the numerical differentiation formulas. Optionally, 
it uses the backward differentiation formulas (also known as Gear’s 
method) that are usually less efficient. It is considered to be good 
for stiff problems while solving a differential-algebraic problem. 
Solver ode23s, a one-step solver is based on an extended Rosen- 
brock formula of order 2. It is could be more efficient than odel5s 
at crude tolerances and is used for the solution of specific type of 
stiff problems where ode!5s are not efficient. 


4. Results and discussion 

4A. Parameter estimation results 

The results of the parameter estimation indicated a good fit be¬ 
tween the model and the measured data. Parameters with low sen¬ 
sitivity and low variability were not further optimized after 
carrying out sensitivity analysis using Aquasim 2.1 d software and 
comparing with previous results of similar cases. Previous study 
indicated that the effect of these parameters on the model out¬ 
comes is quite limited (Batstone et al., 2002a). Optimization of 
the low sensitivity parameters was not carried out since further 
tuning of these parameters requires highly accurate experimental 
data, and changes in the model outputs would be relatively small 
as compared to errors in the data. Those typical parameters were 
selected for the optimization in this work (i.e., similar to the ones 
generally selected for the optimization). 
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Table 3 

Model parameter estimation results with high sensitivity and high variability. 


Parameter 

Name 

Value 

Unit 

^dis 

Dissociation coefficient 

0.75 

d _1 

km-pro 

Monod largest absorption rate of 
propionic acid 

13 

COD/(CODd) 

ks-pro 

Half saturated coefficient of propionic acid 

0.2 

kgCOD/m 3 

/<m-ac 

Monod largest absorption rate of acetic 
acid 

18 

COD/(CODd) 

/Cs-ac 

Half saturated coefficient of acetic acid 

0.2 

kgCOD/m 3 

I<s-h 2 

Half saturated coefficient of H 2 

5 x 10“ 5 

kgCOD/m 3 


In the estimation procedure, the disintegration constant was 
first estimated by matching the model outputs with the measured 
outputs. Then the model outputs for concentrations of acetic 
acid and propionic acid were changed by changing the half satura¬ 
tion constants and maximum uptake rates. Parameters k m , I< s for 
acetic acid and propionic acid were optimized together since they 
exhibit the lowest correlation and the highest relevance (Batstone 
et al., 2003). The optimum values for the parameters are given in 
Table 3. 


4.2. The simulation result of COD 

The raw wastewater COD had a concentration of 15000- 
20000 mg/L, corresponding to an OLR of 30-40 kgCOD/m 3 d. An 
average COD removal rate of 41.7% was achieved by the CSTR, 
while simulated COD removal rates averaged at around 30.0%, 
showing a deviation ratio of 10.1-24.6% between the measured 
and predicted values (Fig. 2). Comparing to the COD removal rates 
(20-30%) in acidogenic phase reported by Li et al. (2006), the pres¬ 
ent CSTR system exhibited a superior performance in COD removal. 
The inconsistency is likely owing to the difference of wastewater 
composition as TCM wastewater was used in the present study. 
The deviation ratios between the measured and predicted values 
indicate that the coefficients and parameters obtained by previous 
work might not be suitable for the acidogenic phase simulation of 
TCM with ADM1 model. 

In the two-phase anaerobic process, wastewater COD is mainly 
removed in the methanogenic phase. The effluent of UASB stabi¬ 
lized at around 1400 mg/L, corresponding to a COD removal rate 
of 86.7%. Only a COD removal efficiency of 65% was achieved in a 
similar, but one-phase UASB reactor while treating chemical syn¬ 
thesis based pharmaceutical wastewater (Ince et al., 2002). Results 



Fig. 2. Measured COD value of CSTR comparison with model simulated value. 
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Fig. 3. Measured COD value of UASB comparison with model simulated value. 
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of this study show that the TPAD system is much superior com¬ 
pared to one-phase processes. It is probably owing to the fact that 
complicated organics were first broken down and converted to 
easily biodegradable components in the CSTR. The experimental 
results did indicate a substantially increase in ethanol concentra¬ 
tion of CSTR effluent. As shown in Fig. 3, the model was able to pre¬ 
dict the COD of UASB effluent with a considerable accuracy. It was 
found that the deviation ratio of the measured and simulated val¬ 
ues fluctuated between -4.8% and 5.0%. 

4.3. The simulation result of VTA 

In the two-phase anaerobic process, the function of acidogenic 
phase reactor is mainly to enhance wastewater biodegradability 
through decomposing large molecular compounds into small 
molecular products. In the present study, VFA and ethanol concen¬ 
trations of CSTR effluent increased substantially, suggesting the 
acidification occurring in the CSTR. It was found an average 
increase in acetic acid from 241 ±1.6 mg/L to 873 ± 2.8 mg/L, in 
propionic acid from 63 ± 1.2 mg/L to 174 ± 1.7 mg/L, in butyric acid 
from 68 ± 0.9 mg/L to 340 ± 2.5 mg/L, valeric acid from null to 
150 ± 2.06 mg/L, and in ethanol from 358 ± 2.1 mg/L to 767 ± 3.2 
mg/L. As a result, the proportions of wastewater VFA and ethanol 
increased from 3.65% to 23.09% after the TCM wastewater reacted 
in the CSTR. Fig. 4 shows the measured and simulated results of 
VFA and ethanol in the CSTR. It can be seen that propionic acid, bu¬ 
tyric acid and valeric acid were well predicted with respective 
deviation ratios of-6.5-7.1%, -2.1-5.7%, and -4.4-4.3%. However, 


acetic acid ethanol was much underpredicted. Compared to the 
measured values, there was a deviation ratio of 13.8-23.2% in the 
simulated results. As a consequence, the VFA content was under¬ 
predicted in this model. It is likely that ethanol present in the 
raw wastewater was converted into acetic acid during the acidifi¬ 
cation process, which was not considered by the model. It was 
noted that Parker (2005) showed that the concentrations of VFAs 
were consistently over-predicted in digesters at short solid 
retention times. These results seems to confirm the previous obser¬ 
vation that it was difficult to accurately predict the VFA concentra¬ 
tion with ADM1 model due to complicated behavior of 
acidification process (Blumensaat and Keller, 2005). 

4.4. The simulation result of pH 

Previous studies indicated that the major products of 
carbohydrates acidogenesis are acetate, butyrate, propionate, and 
ethanol, whose proportion is strongly dependent on culture pH 
(Fang and Liu, 2002). The effluent pH of CSTR was in a range of 
4.8-5.2, and the simulated pH was higher than the measured val¬ 
ues (deviation ratio, -5.8—14.5%), as shown in Fig. 5. It is proba¬ 
bly due to the underpredicted VFA concentrations with the ADM1 
model. 

The effluent pH of UASB fluctuated in a narrow range of 6.5-7.0 
since the acidified intermediates were fermented in the UASB 
(Fig. 6). The model well predicted the effluent pH of UASB with a 
deviation ratio of -2.9-1.4%, indicating that methanogenic bacte¬ 
ria grew well. 



Fig. 4. Measured VAF value of CSTR comparison with model simulated value. 
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Fig. 5. Measured pH of CSTR comparison with model simulated value. 



Fig. 6. Measured pH of UASB comparison with model simulated value. 


5. Conclusions 

ADM1 model was applied to the TCM wastewater treatment 
with TPAD process while considering complete mixing model for 
the CSTR, and axial direction discrete model and mixed series con¬ 
nection model for the UASB. The simulation results indicated that 
the model built was able to well predict the COD removal rate 
(-4.8-5.0%) and pH variation (-2.9-1.4%) of the UASB reactor, 
while failed to simulate the CSTR performance. Comparing to the 
measured results, the simulated acetic acid concentration of the 
CSTR effluent was underpredicted with a deviation ratios of 
13.8-23.2%, resulting in an underprediction of total VFA and COD 
concentrations despite good estimation of propionic acid (devia¬ 
tion ratio, -6.5-7.1%), butyric acid (deviation ratio, -2.1-5.7%) 
and valeric acid (deviation ratio, -4.4-4.3%). It is presumed that 
ethanol present in the raw wastewater was converted into acetic 
acid during the acidification process, which was not considered 
by the model. Additionally, due to the underprediction of acetic 
acid the pH of CSTR effluent was overestimated. 
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